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Abstract 

The generalized Gerasimov-Drell-Hearn (GDH) sum rule is understood based on 
the polarized parton distributions of the proton with the higher twist contributions. 
A simple parameterized formula is proposed to clearly present the contributions of 
different components in proton to Ti(Q 2 ). We find that the contribution of quark 
helicity to T\{Q 2 ) is almost constant (~ 0.123); the twist-4 effect dominates the 
suppression of T^Q 2 ) at Q 2 < 3 GeV 2 , while the twist-6 effect arises a dramatic 
change of Tf (Q 2 ) at Q 2 < 0.3GeU 2 , it implies a possible extended objects with size 
0.2 — 0.3 fm inside the proton. 
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1 Introduction 


The structure of nucleon is an important open question. The perturbative QCD 
(pQCD) at high energy has been proved a successful approach to describe the phenom¬ 
ena in terms of elementary quark and gluon constituents. On the other hand, hadronic 
processes at low transfer momentum are usually described by using meson and baryon 
degrees of freedom. Unfortunately, the transition from quark-gluon to hadron degrees of 
freedom remains shrouded in mystery. 

There is particular interest in the first moment T^Q 2 ) = Jq dxgi(x, Q 2 ) of the spin 
structure functions g\(x, Q 2 ), which has been measured from high Q 2 down to ~ 0 GeV 2 . 
The goal of obtaining universal expressions describing Ti(<5 2 ) at any Q 2 is an attrac¬ 
tive task for both theoretical and phenomenological point of view. In theory, Ti(0) is 
constrained by the Gerasimov-Drell-Hearn (GDH) sum rule [1,2] 

o M 2 k 2 

/ ‘W = *-^- r ‘^ = -T— 0 - 8 - (L1) 

where k is the anomalous magnetic moment of the nucleon. On the other hand, the 
Bjorken sum rule [3] says 


jm [U(Q 2 ) - r;*(<? 2 )] = - 

Q 2 —)■ oo O 


9A 


9v 


( 1 , 2 ) 


this ratio is known very accurately [4]: gA/gv — —1.2695 ± 0.0029. 

The connection of the two sum rules by means of the generalized GDH sum rule is(For 
an overview see, for example Ref. [5]) 


9/l/f 2 

h{Q 2 ) = (1.3) 

which allows us to study the transition between the perturbative partonic structure and 
non-perturbative hadronic picture of nucleon in lepton-nucleon scattering processes. The 
data show that this sum rule at low Q 2 < 1 GeV 2 changes dramatically and exceeds the 
variation bound at higher Q 2 , which has been parameterized (but not explanation) in [6]. 
The explanation of the generalized GDH sum rule is an active subject. For example, the 
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phenomenological constituent quark model [7], the vector meson dominate (VMD) model 
[8], the resonance contributions [9], the chiral perturbation theory (yPT) [10] are used to 
understand the generalized GDH sum rule. 

The purpose of this work is try to expose the partonic structure in the generalized 
GDH sum rule. According to the parton model, the polarized structure functions of 
nucleon at Q 2 > a few GeV 2 and the leading order (LO) approximation are defined by 
the parton distributions as 

g ? GLAP ( X ,Q 2 ) = ^e 2 [5qi(x,Q 2 ) + 8q i (x,Q 2 )], (1.4) 

i 

where the polarized parton distributions obeys the DGLAP evolution equation [11], In our 
previous works [12,13], the DGLAP equation with the parton recombination corrections 
[14,15] are used to derive the unpolarized and polarized parton distributions starting from 
simple initial valence quark distributions at Q 2 ~ 0.064 GeV 2 . Our numerical calculations 
show that the predicted parton distributions consist with the data well both for polarized 
and unpolarized data at Q 2 > 1 GeV 2 . Considering the nonpertnrbative corrections to 
the spin structure function at low Q 2 becoming important, Eq. (1.4) in the full Q 2 range 
should be generalized to 

s?M 2 ) =g? aLAF+ZRS ( x ,Q 2 ) + gr MD (*,Q 2 ) + gr(*,Q 2 ), (1-5) 

where g PGLAP+ZRS (x, Q 2 ) is the solutions of the DGLAP equation with the parton re¬ 
combination (ZRS) corrections at the LL(Q 2 ) approximation [15], while g\ MD (x,Q 2 ) is 
the contributions of the vector meson in virtual photon via the vector meson dominated 
(VMD) model [16], and g RT (x,Q 2 ) is the higher twist corrections. The corresponding 
Erst moment of g\ is 

r?(Q 2 ) = y° glap+zrs (q 2 ) + t\ md (q 2 ) + rf T (Q 2 ). (1.6) 

Since we have known the contributions of g P)GLAP + ZRS an d g\ MD , one can expose the 
properties of Y RT {Q 2 ) after subtracting these two contributions from the experimental 
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data about Fi(Q 2 ). This opens a window to visit higher twist effects at low Q 2 in the 
nucleon structure. We find that (i) the contributions of quark hclicity to Tf(<5 2 ) are 
almost constant (~ 0.123); (ii) the negative twist-4 effect dominates the suppression of 
T^Q 2 ) at Q 2 < 3 GeV 2 , while the twist-6 effect becomes positive at Q 2 ~ 1 GeV 2 and it 
arises a dramatic change of F^Q 2 ) at Q 2 < 0.3Gei/ 2 , these effects can be described by a 
simple parameterized form; (iii) the slope of the higher twist coefficient suddenly enhance 
at Q 2 ~ 1 GeV 2 , which is near an estimation of the constituent quark scale 0.2 ~ 0.3 f m. 

We will give a review of the structure functions in full Q 2 range based on the partonic 
picture in Sec. 2, where the contributions of parton hclicity and the VMD mechanism 
in our previous works are summarized. Then we extract the twist-4 and 6 corrections to 
T^Q 2 ) in Sec. 3, which expose the interesting change of the partonic structure of proton 
at the low Q-scales. The discussions and a summary are presented in Sec. 4. 
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2 Spin structure functions in the full Q 2 range 

We have detailed the first two terms of Eq. (1.5) in our previous two works [13,16]. 
gDGLAP+ZRS are ^ ie con t r ibutions of parton helicities obeying the DGLAP equation with 
the parton recombination corrections [14,15], where the input valence quark distributions 


at p 2 = 0.064GeU 2 have been fixed through a global fitting: 

xu v (x,g 2 ) = 24.3Cte 1 ' 98 (l — x) 2 ' 06 , (2.1) 

xd v (x,g 2 ) = 9.1CU L31 (1 — a;) 3 ’ 80 , (2.2) 

for unpolarized parton distributions, and 

8u v (x, p 2 ) = 40.3x 2 ’ 85 (1 — a:) 2 ' 15 , (2.3) 

5d v (x,p 2 ) = — 18.22a; L41 (l — a;) 4 ' 0 , (2.4) 


for polarized parton distribution. The input distributions Eqs.(2.1-2.4) neglect the contri¬ 
butions of asymmetry sea quark contributions in this work. We consider these corrections 
to Ti are small. 

Using Eqs. (2.1-2.4) we obtain y dglap+zrs (Q 2 ) ~ 0.123, which seems insensitive to 
the complicated behavior of gf at small x [16]. 

Note that 1/p — 0.78 fm consists with a typical proton scale 0.8 ~ 1 f m. We assume 
that [-I is a minimum average transverse momentum of the partons in the proton due to the 
uncertainty principle. Thus, we assume that the Q 2 -dependence of parton distributions 
are freezed at Q 2 < g 2 . Using this assumption we avoid the un-physical singularities at 
Q ~ A qcd- The dashed curve in Fig. 1 is our predicted T PGLAP+ZRS (Q' 2 ). We find that 
it is almost irrelevant to Q 2 . 

The second term of Eq. (1.5) is the the contributions from the interaction of p- and 
cj-mesons in the virtual photon with the proton and they are described by the vector 
meson dominance (VMD) model [17], 
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(2.5) 


xgX MD (x,Q 2 ) = —'£ m ” A ' 7,, ’' W 


8tt V ll(Q 2 + ™l) 2 

In this formula the constants 7 „ are determined by the leptonic widths of the vector 
mesons and m v denotes the mass of p and u at Q 2 < 1 GeV 2 . The cross-sections Aoy(s) 
are the unknown cross sections for the scattering of polarized mesons and nucleons, s is 
the CMS energy squared for the yp collision and s ~ Q 2 jx. In work [16] we assume 


9 1 


VMD, 


(in 2 ')2-A+e / ^)2' 1 A 

_ u y n v) 


x,Q 2 ) = B- 


Q } [(—)~ x+t e(x 0 -x) + 


(Q 2 + rriy ) 2 - x 0 


^ t ^9{x-Xq)\, (2.6) 

In xq 


A = 1 and e is a small positive parameter due to the requirement of integrability of g[ at 
x —)• 0. Thus, we have 


Y MD {Q 2 )= f 1 dxgX MD (x, Q 2 ) ~ A- 
J 0 


m 2 v Q 2 


(D2 1 502 ’ ( 2 - 7 ) 

(■ Q 2 + rn 2 ) 2 

where we take A = 0.162, which corresponds to e = 0.018. The solid curve in Fig. 1 is 
v dglap+zrs^q 2 ) + t vmd^q 2 ^ it shows a wea ker Q 2 -dependence at Q 2 > 1 GeV 2 . 

Comparing the solid curve with Ti(Q 2 ) data. [20-22] in Fig. 1, one can expect that the 
higher twist corrections g^ T play a significant role at low Q 2 to the general GDH sum 
rule. We will detail them in next section. 
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3 Higher twist contributions to the GDH sum rule 

According to the operator product expansion (OPE), the appearance of scaling viola¬ 
tions at low Q 2 is related to the higher twist corrections to moments of structure functions. 
Higher twists are expressed as matrix elements of operators involving nonperturbative in¬ 
teractions between quarks and gluons. The study of higher twist corrections gives us a 
direct insight into the nature of long-range quark-gluon correlations. The higher twist 
corrections to g\ have several representations. The initial parton recombination correc¬ 
tions have been examined in our works [12,13]. The VMD model was used to simulate a 
special higher twist effect in our previous work [16]. In this work, we will try to expose 
the remaining power suppression corrections to Tf. For this sake, we make Tf(Q 2 ) (i.e., 
the data points [28-20] in Fig. l)-[Vp LAP+ZRS (Q 2 ) + Y\ MD (Q 2 )} (i.e., the solid curve in 
Fig. 1). Figure 2 shows such a result at Q 2 > 0.2 GeV 2 , which has been smoothed with 
minimum p j D .o.f.. 

To expose the possible physical information of the curve in Fig. 2, according to QCD 
operator product l/Q 2 -expansion, 


pffT/n2\ _ V'' V2i(Q 2 ) 

1 1 W J-L Q2i-2 ’ 

we take first three approximations 


rf T( 1 ) (Q 2 ) = pfl, 


(3.1) 


(3.2) 


rf T( 4 + 6 ) (Q 2 ) = 


h6 (Q 2 ) + ah (Q 2 )Q 2 
Q 1 


(3.3) 


pHT(4+6+8)^g2^ = h8(Q 2 ) + ^(Q 2 )Q 2 + A4 (Q 2 )Q l ^2 

Q 6 

Then we plot the curves Q 2 rf T( 4 ) (Q 2 ), Q a Y rt(a+q \Q 2 ) and Q 6 rf T( 4 + 6 + 8 ) (Q 2 ) in Fig. 
3. There are following interesting properties of these results: 

(i) <5 6 rf 7(1+6+8 ^((5 2 ) — > 0, if Q 2 —0. This implies that yU 8 vanishes if it is independent 
of Q 2 . Therefore, T( ’ 1+6 \Q 2 ) is an appropriate approximation. 
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(ii) Three curves in Fig. 3 cross at a same point Q 2 ~ 1 GeV 2 . Particularly, the 
intercept jMi of the line, which cuts Q 4 T^ T ^ +6 \Q 2 ) suddenly changes its sing from -0.03 
at Q 2 > 1 GeV 2 to 0.01 at Q 2 < 1 GeV 2 . This result exposes that the correlation among 
partons in the proton has an obvious change near Q ~ 1 GeV. 

(iii) Although Y^ T (Q 2 ) ~ T^ t< ' 4: \Q 2 ) = —0.05 /Q 2 is an acceptable approximation at 
Q 2 > 1.5 GeV 2 , the correct choose is Y^ T (Q 2 ) = Y[ 11{ l+6 \Q 2 ). We use 

T Hn M)(Q2) = _j^_ + _^_ atQ 2 <0 3GeV 2 ( 3 . 5 ) 

to fit the data at Q 2 < 0.3GeP 2 , where we add a parameter e to remove the unnatural 
singularity at Q 2 = 0. The value of e is sensitive to J(0). We find that /x 4 = —0.154, 
)U 6 = 0.037 and e 2 = 0.32GeP 2 . On the other hand, the negative twist-4 effect dominates 
the suppression of Yi(Q 2 ) at Q 2 < 3 GeV 2 , while the twist-6 effect change its sing at 
Q 2 = 1 GeV 2 and its positive effect arises a dramatic change of Yi(Q 2 ) at Q 2 < 0.3 GeV 2 . 

In summary, 


m 2 G 2 

Y P i(Q 2 ) = 0.123 + 0.162 ^ + Tf T (g 2 ) (3.6) 

(■Q 2 + m 

where the HT contributions are 

[ -mmi _ at g2 > 1GeV 2 

Tf T (g 2 ) = J - a08 g/ M2 + °- 01 q 4 /m2 at 0.3 < Q 2 < 1 GeV 2 (3.7) 

0.156/Af- 0.037/M 2 . n 2 p ^ CpV 2 

\ Q 2 +0.32 ” l ” (Q 2 +0.32) 2 11 O 

where M 2 = 1 GeV 2 is the Borel parameter used in the sum rule. We present the com¬ 
parison of onr Tf(g 2 ) with the data [18-20] in Fig. 4. The corresponding Jf(g 2 ) using 
Eqs. (1.3), (3.6) and (3.7) is presented in Fig. 5. 



4 Discussions and summary 


(i) The parton-hadron duality was first noted by Bloom and Gilman [21] in deep inelas¬ 
tic scattering (DIS) and has been confirmed by many measurements. At low energies (or 
intermediate Bjorken variable x and low Q 2 ) DIS reactions are characterized by excitation 
of nucleon resonances; while at high virtuality such processes have a partonic description. 
The smooth high-energy scaling curve essentially reproduces the average of the resonance 
peaks seen at low energies. Burkert and Ioffe [6] indicated that the contribution of the 
isobar A(1232) electroproduction at small Q 2 can describe the general GDH sum rule, 
and they gave 


rL = -0.056 ~ -0.063, at Q 2 = 0.3 - 0.8 GeV 2 (4.1) 

M z 

= 0.010 ~ 0.011, at Q 2 = 0.3 ~ 0.8 GeV 2 , (4.2) 

which are compatible with our prediction Eq. (3.7). 

Our results in Fig.3 indicates the slope of the twist coefficient fi^Q 2 ) in Eq. (3.7) 
suddenly enhance at Q 2 ~ 1 GeV 2 . It implies that the correlation among partons become 
stronger at scale ~ 0.2 fm. We noted that Petronziol, Simula and Ricco [22] reported that 
the inelastic proton data obtained at Jefferson Lab exhibit a possible extended objects 
with size of ~ 0.2 — 0.3 fm inside the proton. 

As we have shown that the twist N > 6 corrections are neglected at low Q 2 . This 
supports our generalized leading order (GLO) approximation [13,16]: if the lower order 
contributions are compatible with the experimental data at low Q 2 , one can image that 
the higher order corrections are cancelable each other. 

In summary, we assume that the parton distributions are still available at Q 2 < 
1 GeV 2 . Using the DGLAP equation with the parton recombination corrections, we in¬ 
dicate that the contributions of Q 2 -dependent parton distributions to the lowest moment 
of the spin-dependent proton structure function are almost constant in full Q 2 range. In 
the meantime, the contributions of a VMD-type nonperturbative part are isolated. After 


9 



removing the above two contributions from the existing experimental data for Tf (Q 2 ), the 
higher twist power corrections present their interesting characters: parton correlations at 
Q 2 ~ 1 GeV 2 and ~ 0.3 GeV 2 show two bend points, where the twist-4 effect dominates 
the suppression of rf(Q 2 ) at Q 2 < 3 GeV 2 , while the twist-6 effect arises a dramatic change 
of r^(Q 2 ) at Q 2 < 0.3 GeV 2 . Within the analytic of these results, we are able to achieve 
a rather good description of the data at all Q 2 values using a simple parameterized form 

of r?(Q 2 ). 
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Figure 1: Contribution of quark hclicity y pglap+zrs (Q 2 ) (dashed curve) and combining 
VMD contribution r? GLAP+ZRS (Q 2 ) + Y\ MD (Q 2 ) (solid curve). Data are from from 
Hermes experiment at DESY [18], the E143 experiment at SLAC [19] and the EGla 
experiment using the CLAS detector at JLab [20]. 
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